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The finite-size scaling method in the equilibrium Monte Carlo(MC) simulations and the finite- 
time scaling method in the nonequilibrium-relaxation simulations are compromised. MC time data 
of various physical quantities are scaled by the MC time data of the dynamic correlation length, 
which corresponds to changing the system size in the finite-size scaling method. This scaling method 
is tested in the three-dimensional ferromagnetic Ising spin model and in the three dimensional ± J 
Ising spin-glass model. The transition temperature and the critical exponents, n and v, are obtained 
by the nonequilibrium relaxation data of the susceptibility and the dynamic correlation length apart 
from the dynamic exponent. We also comment on the definition of the dynamic correlation length 
in the nonequilibrium relaxation process. The Ornstein-Zernike formula is not always appropriate. 



I. INTRODUCTION 

A Monte Carlo(MC) simulation and the finite-size- 
scaling analysis are very important tools in the study 
of phase transitions X The applications serve as a strong 
bridge between the experimental and the theoretical 
physics. Using this method, we may estimate various 
physical parameters, predict unknown properties, and 
propose new experiments on real materials. One of suc- 
cessful applications may be the study on the phase tran- 
sition. We perform the finite-size-scaling analysis on the 
numerical data obtained by the MC simulations. Then, 
the critical temperature and the critical exponents are 
estimated. 

The nonequilibrium relaxation method^— is an alter- 
nate version of the MC simulation that studies the phase 
transitions. The method directly deals with the MC re- 
laxation functions of physical quantities. In this scheme, 
time and size are related by the dynamic scaling relation: 
t ~ ^ Z , which connects the correlation time r and the 
correlation length £ via the dynamic exponent z. In the 
standard equilibrium simulations, we take the infinite- 
time (equilibrium) limit first. Then, infinite-size limit is 
taken by the finite-size-scaling analysis. This procedure 
is reversed in the nonequilibrium relaxation method. We 
take the infinite-size limit first by preparing a very large 
system and stopping the simulation before the finite-size 
effects appear. Then, infinite-time (equilibrium) limit 
is taken by the finite-time-scaling analysis. The critical 
temperature and the critical exponents arc estimated in 
the same manner as the standard finite-size-scaling anal- 
ysis. 

An advantage of the nonequilibrium relaxation method 
is a computational efficiency. We use the nonequilibrium 
MC data that are usually discarded in the equilibrium 
simulations. It is very important particularly when the 
computational resources at hand arc limited. 

One shortcoming in the nonequilibrium-relaxation 
method is an ambiguity of the dynamic exponent, z. 
We obtain the transition temperature and the critical 
exponent, 7, by the finite-time-scaling analysis on the 
susceptibility^— However, the exponent v is not solely 



obtained but it appears as a form zv. We need to esti- 
mate z independently in order to obtain v. The dynamic 
exponent takes a large value in the frustrated and/or ran- 
dom systems. It sometimes depends on the temperature 
in spin-glass modelsJ^ The numerical estimate of z is a 
tough task in these simulations. 

In this paper, we introduce another scaling method 
based on the nonequilibrium relaxation scheme. The 
method is free from the estimate of z. We focus on 
the nonequilibrium relaxation function of the correlation 
length, which is called the dynamic correlation length. It 
serves as the system size in the finite-size-scaling analy- 
sis. The nonequilibrium relaxation functions of various 
physical quantities are scaled by the dynamic correlation 
length. Wc obtain the transition temperature and the 
critical exponents by this scaling analysis in the same 
manner as the finite-size scaling. We just replace the 
system size L of the finite-size scaling analysis with the 
dynamic correlation length £ (t) , where t is the MC time 
step. We also propose a definition of the dynamic cor- 
relation length. The Ornstein-Zernike formula^ is often 
used to define this value. We show that this definition 
is not always appropriate particularly in the nonequilib- 
rium relaxation regime. 

Section |TT] explains models and a definition of the dy- 
namic correlation length. Section Hill explains the scaling 
procedure. Numerical results are presented in Sec. IIVI 
The ferromagnetic Ising model and the ±J Ising spin- 
glass model are taken as examples. Section [V] is devoted 
to the summary and discussions. 

II. MODEL AND DEFINITION OF THE 
DYNAMIC CORRELATION LENGTH 

Let us consider the following spin model to demon- 
strate the present scaling method. 

W = — ^ JijSiSj, 

where denotes the nearest-neighbor pairs in the cu- 
bic lattice, Jij denotes the exchange interaction, and Si 
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denote the Ising spin. We consider the uniform ferromag- 
netic model (Jij = 1) and the spin- glass model (Jy = ± J 
with an equal probability) in this paper. 

The correlation length of an ordered domain in spin 
models is usually defined by the correlation function 
C(r) = (SiS i+r ) asH 



C{r) = (SiS. l+r ) oc 



exp[-r/£] 



in the equilibrium state near the transition temperature. 
In the Monte Carlo simulations, it is easier to estimate 
the correlation length by the Fourier transform of the 
susceptibility, x(k), as ~ 
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Here, fc m in denotes the smallest wave vector along one 
direction in a finite-size lattice. The expression is exact 
in the limit of k m - m — > 0. 

The 1/r-correction form in the spin correlation func- 
tion is only guaranteed by the mean-field theory in the 
equilibrium stated It is not trivial that we may use 
Eq. ([TJ in the nonequilibrium relaxation process. The 
system size L is much larger than the correlation length 
in the nonequilibrium process. In such case, the spin 
correlation function may decay in a simple exponential 
form without the 1 /r-corrcction. In order to check the in- 
dependence, we plot the correlation functions, C(r) and 
rC(r), against r in the nonequilibrium regime (t = 501) 
and in the equilibrium regime (t = 3981). The system is 
the uniform ferromagnetic model with Jy = 1, the sys- 
tem size L = 299, and the temperature T = 4.515, which 
is above the transition temperature, T c = 4.511525^ 
Figure Q] shows the result. The linearity of rC(r) is bet- 
ter than that of C(r) when t = 3981. The 1/r-correction 
is necessary in estimating the correlation length in the 
equilibrium state. On the other hand, the linearity is 
poor and convex upward when t = 501. The linearity 
of C(r) seems to be better. If the correlation function 
exhibits the simple exponential decay as 

C(r) = {SiSi+ r )<xexp[-r/£\, 
the correlation length should be estimated by 
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It is difficult to judge a use of Eq. ([T]) or Eq. ^ just by 
the linearity of C (r) or rC (r) . The difference is small as 
shown in Fig. [TJ Therefore, we compare the time depen- 
dences of the dynamic correlation length estimated by 
both definitions. We judge two definitions by the scaling 
behavior. 

Figures [5] and [3] compare the estimates of using 
both definitions, Eqs. (JlJ and ©. The dynamic corre- 
lation length of the size-L system, L) divided by the 
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FIG. 1: (Color online) r dependence of the spin correlation 
function in the nonequilibrium regime (t = 501) and in the 
equilibrium regime (t = 3981). The system is the uniform 
ferromagnetic model in three dimension. The linear system 
size L = 299, and the temperature T = 4.515. The lower two 
lines are plots of C(r), and the upper two lines are those of 
rC(r). 




FIG. 2: (Color online) The dynamic correlation length 
£(i, L) over L is plotted against t 1//z over L. Data labelled 
by (1) are estimated by Eq. (1), and those labelled by (2) 
are estimated by Eq. (2). The system is the uniform 
ferromagnetic Ising model in three dimensions. The 
temperature T = 4.512 is very close to the transition 
temperature. The estimate of the dynamic exponent 
2 = 2.05 is taken from the Ref.— . 



linear scale L is plotted against t x l z /L. We expect that 
t;(t, L) oc t x l z from the dynamic scaling relation. The 
correlation length divided by L becomes scale invariant 
at the transition temperature. Therefore, the present 
plot should ride on the single line with its slope unity, if 
we choose a correct value of z that satisfies the dynamic 
scaling relation. 

In the uniform ferromagnetic Ising model (Fig. [2]), the 
dynamic correlation length estimated by Eq. ^ yields a 
good scaling and the slope unity. The dynamic exponent 
z = 2.05 is taken from the estimate given by Ito et al^. 
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FIG. 3: (Color online) The dynamic correlation length 
£(t,L) over L is plotted against i 1//z over L in the ±J Ising 
spin-glass model in three dimensions. Data labelled by (1) 

are estimated by Eq. (1), and those labelled by (2) are 
estimated by Eq. (2). Three temperatures, T — 1.4, 1.18, 
and 1.10, are above, near, and below the transition 
temperature, respectively. We estimated the dynamic 
exponent so that the dynamic correlation length are 
well-scaled at each temperature. 



The finite-size effects appear as the system approaches 
the equilibrium state. 

The dynamic correlation length given by Eq. ([TJ is 
well-scaled but the slope does not agree with unity. If 
we choose z so that it yields the slope unity, the relax- 
ation data of different sizes do not ride on the same line. 
The slope suggests that £(*, L)/L oc (t x / z /L) 1 - 13 , but the 
scaling fails if we take z = 2.05/1.13 = 1.81 for t 1 / z /L. 
Therefore, the dynamic scaling relation contradicts itself 
even in the ferromagnetic model. It is an evidence that 
a use of the Ornstcin-Zcrnikc expression, Eq. ([T]), is not 
appropriate in the nonequilibrium relaxation scheme. 

Figure [3] shows the same £(t,L)/L scaling in the ±J 
Ising spin-glass model in three dimensions. Relaxation 
data at three temperatures are plotted in the same fig- 
ure by shifting arbitrary in the vertical direction. The 
temperatures are above (T = 1.40), near (T = 1.18), 
and below (T = 1.10) the transition temperature. The 
dynamic exponent at each temperature is estimated so 
that the relaxation data fall onto the same line with the 
slope unity. The linearity is better if we use Eq. © . The 
estimate of £(t,L) using Eq. (|T|) bends upward as the 
system approaches the equilibrium limit, while the slope 
approaches unity in the nonequilibrium limit (t — > 0). 
This bending behavior may mislead us to underestimate 
the value of z. It is also noted that the value of z de- 
pends on the temperature as zT ~ 7.— The finite-size 
effects appear when the temperature is above the transi- 
tion temperature (T = 1.40). The size effects are weak 
in the spin-glass phase (T = 1.18 and 1.10). 

It is shown that the dynamic correlation length in the 
nonequilibrium relaxation process should be estimated by 
using Eq. ([3J . The dynamic scaling relation is satisfied in 



the nonequilibrium process if we take this definition. We 
perform the dynamic correlation-length scaling analysis 
using Eq. ^ in this paper. 

III. SCALING PROCEDURE 

We explain our scaling procedures step by step in this 
section. The first two steps below are the standard pro- 
cedure of the nonequilibrium relaxation method. The 
relaxation functions of the physical quantity and the cor- 
relation length are obtained. 

1. We prepare a system with a linear scale L, and per- 
form a MC simulation at the temperature T. We 
start the simulation from the paramagnetic state. 
The correlation length is zero in the initial state. A 
physical quantity, A, is measured at each MC step, 
t, and is stored in memory as A(t, L, T). Changing 
a random number sequence, we perform indepen- 
dent simulations and take an average of A(t,L,T) 
over these runs. 

2. We change the system size and perform simula- 
tions. Relaxation functions A(t, L, T) of different 
sizes are compared to check the finite-size effects. 
Only data that are free from the finite-size effects 
are used in the scaling analysis. Examples arc Figs. 
4 and 7 in Sec. IIV1 The data free from the size 
effects are denoted by A(t,T). We obtain a relax- 
ation function of the dynamic correlation length, 
£(t,T), as explained in Sec. HT1 

We will perform three scaling analyses below. The 
first scaling analysis determines the dynamic exponent, 
z. Then, the second one determines the anomalous ex- 
ponent r], and the last scaling analysis determines the 
critical temperature, T c , and the exponent v. It is noted 
that the first scaling analysis is not necessary and is inde- 
pendent from the other two scaling analyses. It should be 
done only when we need a value of z. This is a clear differ- 
ence from (and is probably an advantage over) the stan- 
dard nonequilibrium relaxation method with the finitc- 
time-scaling analysis. 

3. The relaxation function of the dynamic correlation 
length should satisfy the dynamic scaling relation: 
£(t,T) ~ t x / z . At the temperature near the crit- 
ical temperature, We plot £(t,L)/L versus t x / z /L 
so that all the scaled data fall on a single scaling 
function. This scaling is possible because £/L is 
scale-invariant at the critical temperature. Only a 
value of z is a control parameter. Examples are 
Figs. 2 and 3. 

Let us take the magnetic susceptibility, Xi as a physical 
quantity for simplicity. 

4. Using the relaxation function of x(t, T) and T), 
we perform the scaling analysis. Since the simu- 
lation starts from the paramagnetic state, both \ 
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and £ are expected to increase algebraically in time 
as x(t,T) ~ £(i,T) 2 ~ r ' when T is near T c . Since 
X oc (, 2 ~ v in the critical region, we expect that the 
nonequilibrium relaxation functions of %(t, L) and 
£,(t,L) depend on L as 
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We plot x{t,L)/L 2 -^ against £(t,L)/L in a log- 
log scale and determine 77 so that all the relaxation 
functions falls on the same line with a slope 2 — 77. 
Examples are Figs. 5 and 8 in Sec. IIV1 It is a 
one-parameter scaling analysis. The determination 
of 77 is easily done first by this scaling plot. Here, 
it is not necessary to know the critical tempera- 
ture precisely. We only need the rough estimate. 
Since we work with the nonequilibrium relaxation 
method, the critical relaxation is observed near the 
critical temperature in the nonequilibrium process. 
It is sufficient to perform this scaling analysis at 
one temperature near T c . 

5. The critical temperature, T c , and the critical expo- 
nent, is, arc determined in the last scaling anal- 
ysis. The equilibrium value of the correlation 
length diverges as £ ~ \T — T c \~ v . At the same 
time the susceptibility x(i, T) should be scaled by 
T) 2 ~ n with 77 estimated above. Therefore, we 
plot the relaxation functions x(t, T)/£(t, T) 2 ~ v ver- 
sus £(t,T)/\T — T c \~ v for various temperatures and 
determine T c and v so that all the data ride on a 
single scaling function. Examples are Figs. 6 and 
9. in Sec. HVl 



It is noted that the finite size L in the finite-size- 
scaling analysis is replaced by the dynamic correlation 
length £(£, T) in the present analysis. We consider that 
the present scaling analysis is a natural extension from 
the finite-size-scaling analysis when we work within the 
nonequilibrium relaxation scheme. 



IV. NUMERICAL RESULTS 



icr 



1(T 



i(r 



i. 10 2 



10 J 



10 u 



L=299 - 


3^ — 






199 ^ 




- 159 




79 






__L- — 




7=4.512 



10 1 



10 z 



l(f 



10" 



FIG. 4: (Color online) Size dependences of raw relaxation 

data of the dynamic correlation length, £(t, L), and the 
magnetic susceptibility, L) in the ferromagnetic model. 
The temperature, T = 4.512, is a little above the transition 
temperature. 
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FIG. 5: (Color online) The \ ~t scaling plot of Eq. © in 
the ferromagnetic model. The raw relaxation data of Fig. [4] 
are scaled. 



1.85 and 1.95, we search for values of T c and v so that 
the scaled data yield the best scaling behavior. The ob- 
tained values are listed in Fig. [5] and are consistent with 
the previous estimates.— 



A. Ferromagnetic model 



B. Spin-glass model 



We present the numerical results in the ferromagnetic 
Ising model in three dimension. Figure @] shows the raw 
relaxation data of the dynamic correlation length and the 
magnetic susceptibility. The dynamic correlation length 
is estimated using Eq. ([2]). We plot in Fig. [5] the scaling 
result that determines 77. The relaxation data are scaled 
with 2 — r\ = 1.90 ± 0.05, where the logarithmic slope of 
the scaled function agrees with 2 — 77. 

The dynamic correlation-length scaling of the suscep- 
tibility is shown in Fig. [B] Here, only the nonequilibrium- 
rclaxation data that are free from the finite-size effects 
are plotted. For a estimated value of (2 — 77) between 



In this subsection, we present the numerical results 
of the ±J Ising spin-glass model in three dimensions. 
The spin-glass correlation length is estimated from the 
Fourier transform of the spin-glass susceptibility, ;^sg- It 
is defined by 
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where the bracket [• • • 
age, and the bracket 



denotes the configurational aver- 
■ ) denotes the thermal average. 
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FIG. 6: (Color online) The dynamic correlation- length 
scaling plot of the magnetic susceptibility in the 
ferromagnetic model. 

The thermal average is estimated by the average over 
different real replicas: 

A=l 

A replica number is denoted by m, and the superscript 
A denotes a replica index. We prepare m real replicas 
for each random bond realization with different initial 
spin configurations. Spin states of each replica are up- 
dated using different random number sequences. The 
thermal average is taken only by this replica average in 
our nonequilibrium relaxation scheme. A replica number 
controls an accuracy of the thermal average. It is set 
to 256 in this paper. The dynamic spin-glass correlation 
length in this paper is estimated by Eq. ([2]). where x is 
replaced by Xsg- It is usually defined by Eq. ([T]) . 17 '' 18 

Figurc[7]shows the raw relaxation functions of the spin- 
glass susceptibility and the dynamic spin-glass correla- 
tion length. Figure [5] shows the scaling plot determining 
a value of rj. The straight-line scaling is possible near 
and above the transition temperature. The scaling poor 
in the low-temperature phase. Using the estimated value 
of 2 — rj wc plot the dynamic correlation- length scaling re- 
sults in Fig. [5] The scaling is good when the temperature 
is above the estimated spin-glass transition temperature, 
T sg = 1.18. The scaling behavior becomes poor below 
T sg . The estimated critical exponents may be invalid in 
the low-temperature phase. 

The previous estimates^— for T sg lie between 1.1 
and 1.2. They are roughly categorized into two groups. 
On o 18 ' 21 gives T sg close to 1.1 and v close to 2. The 
otheriSi 2 ^ gives T sg close to 1.2 and u close to 1.3. The 
present result is consistent with the latter group. The 
latter group mostly takes the dynamic approach to the 
phase transition. Recently, Hukushima and Campbell 2 ^ 
discussed that this discrepancy can be understood by the 
strong corrections to scaling. 

Campbell et al£^ proposed the /3-scaling analysis, 



FIG. 7: (Color online) Size dependences of raw relaxation 
data of the dynamic correlation length, £ a g(i, L), and the 
spin-glass susceptibility, Xsg(t,L) in the spin-glass model. 

Three temperatures are above, near, and below the 
transition temperature, respectively. Data of T = 1.10 are 
divided by 2, and those of T = 1.25 are multiplied by 2 in 
order to separate from data of T = 1.18. 
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FIG. 8: (Color online) The Xsg-£sg scaling of Eq. ([3j) in the 
spin-glass model. The raw relaxation data of Fig. [7] are 
scaled. A straight line is a guide for eyes with a slope 
2 — 77 = 2.2. 
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FIG. 9: (Color online) The dynamic correlation-length 
scaling plot of the spin-glass susceptibility. 
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FIG. 10: (Color online) The dynamic correlation-length 
scaling plot of the spin-glass susceptibility using \0% — /3 2 | as 
the scaling variable. 



which uses (/3 2 — /3c) as the scaling variable. They es- 
timated the transition temperature and the critical ex- 
ponents as T c = 1.11, v = 2.72, and 2 - 1) = 2.40(4). We 
also try this ^-scaling analysis in Fig. [TUl Our estimates 
are T c = 1.11, v = 2.62, and 2 — r\ = 2.35, which are con- 
sistent with their estimates but disagree with our present 
estimates using (T — T c ). This discrepancy suggests that 
the present numerical simulation is not sufficient to ex- 
tract the true critical phenomena both in a size scale 
and a time scale in the spin-glass model^ We checked 
that the /3-scaling is possible using the same transition 
temperature and critical exponents in the ferromagnetic 
model. (Figures are not shown.) The present size and 
time scales are sufficient in the ferromagnetic model. 



the scaling analysis. The examples are Figs. [5] and [H Al- 
though the raw relaxation data (Figs. 2] and [7]) exhibit 
the finite-size effects, we can scale them to one scaling 
line without the size effects. The critical divergence is 
observed from very early stage of the nonequilibrium re- 
laxation process, if we scale the data by the dynamic 
correlation length. 

The present dynamic-correlation-length-scaling analy- 
sis is regarded as an extension of the finite-size-scaling 
analysis replacing the size L with the dynamic correla- 
tion length £(t). We may consider that the finite-time 
relaxation data at a time t corresponds to the equilib- 
rium data of the size L with L = £(i). 

Since the present scaling method is entirely based on 
the dynamic correlation length, the definition is very im- 
portant. We found that a use of the definition based on 
the Ornstein-Zcrnikc formula is not appropriate in our 
scheme. Since the formula is based on the mean-field ap- 
proximation, the present nonequilibrium process may be 
out of the applicable range of the approximation. The re- 
laxation function of the dynamic correlation length esti- 
mated by Eq. ([T]) exhibits an extra increase before reach- 
ing an equilibrium value, as shown in Figsj2]and[3] The 
logarithmic slope, which is 1/z, is then overestimated. It 
affects a value of v in the conventional finite-time-scaling 
analysis of the nonequilibrium relaxation method^ The 
introduced definition, Eq. ([2J , is based on a simple expo- 
nential decay of the correlation function. The relaxation 
function exhibits a normal behavior. It algebraically in- 
creases with t x / z , and eventually converges to a finite-size 
value. We also comment that the definition Eq. ([2]) may 
be used even in the equilibrium state where the mean- 
field approximation is not valid. 



SUMMARY AND DISCUSSION 
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